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*6.  Abstract 

The  normal  venting  capacity  of  marine  cargo  tanks  appears  to  be  adequate  for 
relieving  vapor  generated  when  Integral  tanks  are  exposed  to  an  external  fire. 
However,  the  unwetted  tank  «mlls  surrounding  the  ullage  are  subject  to  sub- 
stantial weakening  In  those  cases  where  the  fire  heat  flux  increases  the  wall 
temperature  above  lOOOtp.  Convection  and  radiation  can  adequately  cool  the 
unwetted  wall  only  for  incident  heat  flux  less  than  16,000  BTD/hr/ft^.  Thus, 
the  failure  mode  of  cargo  tanks  with  an  external  fire  may  be  due  not  to  the 
vent  system,  but  rather  to  heating  of  the  unwetted  wall. 
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1.  Pg.  4,  Par.  1,  sentence  6:  Replace  20  pslg  with  6. 


2.  Pg.  11,  under  section  4:  Delete  third  and  fourth  sentences  of 
the  first  paragraph.  Add  the  following  to  the  first  paragraph: 

Using  a modified  theory  of  Clarkson  (1956) , estimates  were  made 
of  the  minimum  internal  pressure  loadings  required  to  initiate  failure 
of  cargo  tanks  of  three  representative  vessel  designs.  For  the  offshore 
barge,  an  estimate  of  6 pslg  was  obtained;  for  the  Inland  barge,  4 psig: 
and  for  a large  tankshlp,  8 psig.  The  average  value  for  these  three, 

6 pslg,  will  be  utilized  in  this  report  as  a nominal  failure  pressure  for 
tank  vessel. 

3.  Pg.  13,  fifth  line:  replace  "20  pslg"  with  "6  pslg". 

1 sixth  line:  replace  "20  psig  or  3.2  psig"  with  "6  psig  or  1 psig". 

i 4.  Pg.  27,  third  line  from  bottom:  replace  "20  pslg"  with  "6  psig". 
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6.  Pg.  41,  last  sentence:  replace  existing  sentence  with  the  following: 
"In  general,  our  results  suggest  that  the  exposed  wall  area  should  not 
be  more  than  1000  times  vent  area  if  the  pressure  is  not  to  exceed  6 pslg 
[ for  a typical  cargo  and  heat  flux". 
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SUMMARY 


RECCMIMEMDATIOMS 


Venting  systems  are  Installed  on  marine  cargo  tanks  to  relieve  pressure 
differentials  which  arise  between  the  vapor  space  and  the  ambient.  The 
capacities  of  the  venting  systems  are  designed  to  accommodate  the  flow  of 
vapor  and  air  during  cargo  transfer  operations  under  a normal  back  pressure 
of  about  2 psig.  However,  in  the  event  of  a fire  external  to  the  tank,  sub- 
stantial pressure  rise  may  occur  due  to  vapor  generation  and  due  to  heating 
of  the  trapped  vapor.  The  hazard  of  pressure  rise  due  to  external  fire  and 
the  design  requirements  of  cargo  venting  systems  to  control  the  hazard  are 
being  reviewed  by  Arthur  D.  Little,  Inc.  as  part  of  Contract  No.  DOT-CG- 
42,357-A.  Three  other  hazardous  situations  which  occur  in  marine  cargo 
transport  have  been  analyzed  in  earlier  reports:  Flame  propagation  into 
the  vent  system,  blockage  of  venting  relief  valves,  and  overpressure  due 
to  excessive  loading  rates. 

In  the  present  report,  the  external  fire  problem  is  analyzed  taking  into 
account  variations  in  the  following  parameters: 

• Levals  of  heat  flux  incident  on  the  tank; 

• Exposed  areas  (wetted  and  unwetted) ; 

• Tank  wall  thickness  and  thermal  and  structural  properties; 

• Cargo  heat  of  vaporization; 

• Molecular  weight  of  the  cargo  vapor;  and 

• Capacity  of  the  venting  system. 

The  case  of  liquid  venting  (which  might  occur  if  the  vessel  was  listing) 
is  not  Included  in  this  analysis,  but  has  been  studied  for  high  venting  rates 
in  an  earlier  report. 

Numerical  methods  were  developed  for  predicting  the  pressure  rise  and 
unwetted  wall  temperature  as  a function  of  time  for  any  apecifled  combina- 
tion of  parameters  listed  above.  By  exercising  this  powerful  numerical 
tool,  it  was  possible  to  resolve  such  issues  as: 

(1)  Does  a tank  exposed  to  fire  fail  due  to  Inadequate  venting 
capacity  or  weakening  of  the  heated  walls? 

(2)  What  la  the  longest  fire  duration  which  can  be  tolerated 
before  failure? 

(3)  What  venting  capacity  is  required  for  a given  tank/ cargo/ fire 
situation? 
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B.  SUMMARY  OF  FINDINGS 


The  normal  venting  capacity  of  marine  cargo  tanka  appears  to  be  adequate 

for  relieving  vapor  generated  when  Integral  tanks  are  exposed  to  an  external 

fire.  However,  the  unwetted  tank  walls  surrounding  the  ullage  are  subject 

to  substantial  weakening  In  those  cases  where  the  fire  heat  flux  Increases 

the  wall  temperature  above  1,000F.  Convection  and  radiation  can  adequately 

cool  the  unwetted  wall  only  for  Incident  heat  flux  less  than  16,000  Btu/ 

2 

hr/ft  . Thus,  the  failure  mode  of  cargo  tanks  with  an  external  fire  may  be 
due  not  to  the  vent  system,  but  rather  to  heating  of  the  unwetted  wall. 

Cargo  boiling  behind  the  wetted  wall  Is  not  the  only  cause  of  pressure 
rise;  heating  of  the  confined  vapor  space  by  the  unwetted  wall  also  Increases 
the  tank  pressure.  These  two  mechanisms  produce  a two  stage  pressure  history 
leveling  off  after  about  20  minutes  from  the  start  of  the  fire.  The  pre- 
dicted pressure  rise  under  representative  flr3  conditions  Is  less  than  2 pslg; 
that  Is,  the  normal  venting  capacity  appears  to  be  adequate  to  release  vapor 
rapidly  generated  by  an  external  fire.  The  following  unusual  conditions 
would  cause  the  tank  pressure  to  exceed  10  pslg  during  a fire: 

• Ratio  of  vent  area  to  the  wetted  wall  area  exposed  to  fire:  Less 
than  .0005. 

2 

• Ratio  of  heat  flux  (Btu/hr/ft  ) to  heat  of  vaporization  (Btu/lb): 
Greater  than  660. 

The  pressure  rise  Is  also  aggravated  by  a large  vapor  space  and  by 
cargo  vapors  of  low  molecular  weight.  Those  combinations  of  variables 
leading  to  pressures  exceeding  4 pslg  can  be  Identified  by  expression  (24) 
which  appears  on  page  46. 

Since  the  analyses  employed  In  this  work  utilize  well  established  for- 
mulas for  venting  and  heat  transfer  and  take  Into  account  all  of  the  Im- 
portant effects  that  may  occur  during  an  external  fire.  It  Is  expected  that 
predictions  of  ullage  pressure  will  Indeed  correspond  closely  to  the  pressures 
encountered  In  actual  fire  situations. 

C.  RECOMMENDATIONS 

Because  the  weakening  of  cargo  tank  walls  by  an  external  fire  can  cause 
catastrophic  failure  of  the  tank,  additional  work  needs  to  be  carried  out  to 
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better  determine  the  precise  conditions  under  which  barge  and  tanker  vessels 
are  excessively  heated  and  to  establish  means  of  prevention.  It  is  recommended 
that  studies  be  conducted  to: 

• Determine  the  duration  and  size  of  marine  fires  resulting  from 
spills,  engine  room  fires,  ignited  vent  gases,  and  other  sources; 

• Establish  the  effect  of  the  heat  flux  of  each  type  of  flame  on  the 
temperature  of  cargo  tank  walls,  taking  Into  account  conduction  to 
the  supporting  structure  and  unexposed  walls;  and 

• Evaluate  the  potential  for  external  Insulation  and  the  effects  of 
double  hulls. 

Such  studies  should  be  accompanied  by  experiments  to  test  the  validity 
of  the  analyses  and  the  feasibility  of  preventive  methods. 


II.  INTRODUCTION 


A.  PROBLEM  STATEMENT 

Normally,  the  air  displaced  from  a tank  during  loading  Is  vented,  along 
with  the  vapor  which  Is  generated  from  the  liquid-cargo  surface.  The  vent- 
ing system  usually  Includes  a relief  valve  and  represents  a constriction 
which  at  normal  venting  rates  causes  the  tank  pressure  to  rise  to  1 - 3 pslg. 
Although  failure  points  vary  widely,  a typical  tank  structure  Is  capable  of 
withstanding  about  f)pslg*(lf  constructed  In  accordance  with  ABS  rules)  and 
Is  therefore  capable  of  handling  the  normal  back  pressure  of  up  to  3 pslg. 

If  there  Is  an  external  heating  source  such  as  a fire,  the  liquid  evapora- 
tion rate  Increases  and  the  back  pressure  Increases  accordingly.  Fire  may 
result  from  collision  with  another  vessel,  resulting  In  the  release  and  Ig- 
nition of  flammable  cargo.  The  tank  walls  heat  up,  the  liquid  cargo  heats 
up  and  evaporates  at  a faster  rate,  and  under  some  conditions,  boiling  may 
occur.  This  rapid  release  of  vapor  with  a finite  ve'^tlng  capacity  causes  a 
continuous  pressure  rise.  To  compound  the  problem,  the  vapor  Is  heated  by 
the  unwetted  hot  tank  walls  and  this  Increases  the  tank  pressure  even  more. 
Under  extreme  conditions  the  back  pressure  may  Increase  to  the  tank  failure 
point  (which  Is  lower  than  normal  because  of  external  heating).  The  purpose 
of  our  analysis  Is  to  determine  the  conditions  under  which  tank  failure  could 
occur,  and  to  estimate  the  elapsed  time  before  failure.  Key  variables  are 
the  vent  area,  the  exposed  wall  area,  the  fire  characteristics,  and  the  li- 
quid cargo  properties. 

Our  approach  has  been  to  perform  a detailed  transient  heat  transfer 
analysis  of  a partlally-full  rectangular  tank  exposed  to  lateral  heating. 

The  vapor  generation  rate  and  transient  wall  temperature  obtained  from  this 
analysis  are  then  used  to  describe  the  pressure  rise  and  decreasing  wall 
strength,  respectively.  As  shown  In  Figure  Il-l,  the  cross-over  of  the  In- 
creasing pressure  and  decreasing  strength  Indicates  tank  wall  failure. 

An  earlier  analysis  by  the  National  Research  Council  (1973)  contained 
a comprehensive  survey  of  the  heat  flux  to  be  expected  from  fires.  The  N.R.C. 
report  also  derived  relief  valve  venting  capacities  by  assuming  that  in  steady 

state  all  of  the  heat  flux  goes  to  generating  vapor.  However,  several 

_ 

Based  on  plastic  deformation  theory  applied  to  actual  bulkhead  thickness 
and  stiffener  spaclngs  (see  Vent  System  and  Loading  Criteria  for  Avoiding 
Tank  Overpressurlsatlon,  part  3 of  Contract  DOT-CG-42,357-A). 
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• Boiling  of  Cargo  at  Wetted  Wail 
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FIGURE  11-1  METHOD  OF  ANALYSIS 
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limitations  make  the  N.R.C.  study  Inadequate  to  describe  the  transient  re- 
sponse of  an  Integral  tank  to  external  fire:  (1)  the  "critical  flow"  venting 
assumption  (Equation  (28),  p.  51)  is  not  appropriate  for  low  pressure  venting; 
(2)  the  N.R.C.  study  did  not  Include  transient  heatup  of  the  tank  wall;  (3) 
secondary  heating  of  the  vapor  space  was  not  included;  and  (4)  specific 
values  of  liquid  height  above  sea  level  and  unwetted  wall  height  were  not 
considered;  Instead,  general  factors  were  employed  to  represent  the  tank 
wall  area,  exposure  factor,  and  atmospheric  attenuation  factors.  In  these 
respects  the  present  analysis  is  more  general  and  more  readily  applicable 
to  Integral  tanks  in  barges  and  chemical  tankers. 


B.  FIRE  SITUATIONS 


Four  fire  situations  were  considered  In  formulating  the  problem  (see 
Figure  II-2) ; 

• Fire  In  the  hull; 

• Fire  on  deck; 

• Adjacent  fire;  and 

• Spot  fire. 

The  “adjacent"  fire  was  selected  as  the  most  useful  case  to  analyze  because 
(a)  the  fire  duration  Is  great  compared  to  the  spot  or  deck  fire,  (b)  the 
frequency  of  occurrence  Is  relatively  great,  (c)  the  fire  size  Is  larger 
than  the  spot  fire  and  large  enough  to  envelop  the  entire  side  wall,  and 
(d)  both  wetted  and  unwetted  portions  of  the  tank  wall  are  heated.  A review 
of  DOT  casualty  reports  covering  a six  year  period  showed  that  the  "adjacent" 
fire  and  the  engine  room  fire  were  comparatively  more  often  encountered 
than  other  types,  with  several  casualties  cited  annually. 

The  hull  fire  requires  adequate  ventilation  below  the  tank  and  occurs 
with  Independent  tank  construction  rather  than  with  Integral  tanks.  This 
situation  has  been  analyzed  by  AOL  In  a recent  study  of  pressurized  liquid 
propylene  tanks  (Raj  [1975]).  The  deck  fire.  If  prolonged.  Is  more  likely 
to  cause  simple  melting  of  the  tank  wall  than  fracture  due  to  Inadequate 
venting,  since  the  vapor  generation  rate  (which  absorbs  heat  and  protects 
the  wall)  Is  relatively  low. 

il 
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Fire  in  Huil 

Fire  on  Deck 

Adjacent  Fire 

Spot  Fire 

Adjacent  tank  punctures 
or  fails,  spilled  liquid 
ignites. 

Flammable  cargo,  cargo 
vapors,  or  other  liquid  ig- 
nited by  spark  or  smoking 
material. 

Adjacent  vessel  suffers 
cargo  spill  or  explosion 
or  cargo  spill  ignites. 

Fire  in  engine  room  trans- 
mitted through  seal  or 
bulkhead. 

SuHace  exposed  to  heat  flux 


C.  RANGES  OF  KEY  PARAMETERS 


1.  Heat  Flux  from  External  Fire 

The  Incident  flux  from  hydrocarbon/air  fires  has  been  measured  and  cal- 
culated by  a number  of  Investigators  for  widely  varying  conditions.  Values 

2 

range  up  to  90,000  Btu/hr/ft  , and  depend  on  the  optimal  beam  length, 
carbon-hydrogen  ratio  of  the  fuel,  and  separation  distance.  The  USCG  regu- 
lation (46  CFR  54  [1971])  assumes  a value  of  34,500  Btu/hr/ft^,  which  is 
reasonable  since  it  is  comparable  to  the  heat  flux  to  the  furnace  walls  of 
an  oil-fired  utility  boiler.  The  National  Research  Council  (1973),  p.  24  ff, 

compiled  data  for  a number  of  fuels  and  fire  sizes,  supporting  a representa- 

2 2 
tive  value  of  34,500  Btu/hr/ft  with  a range  from  5,400  to  47,540  Btu/hr/ft  . 

2 

Hot  spots  have  been  observed  receiving  up  to  90,000  Btu/hr/ft  , according  to 
Anderson  and  Streslno  (1963).  Based  on  this  data  base,  we  employed  the  fol- 
lowing values  of  heat  flux: 

Heat  Flux  Values 

2 

Representative  Value  32,000  Btu/hr/ft 

Range  of  Values  16,000  to  64,000  Btu/hr/ft^ 

Recent  studies  by  Markstein  (1974)  and  Lee  (1975)  on  small  flames  confirm 
that  the  major  heat  flux  component  is  radiative. 

2.  Tank  and  Vent  Configuration 

For  the  purposes  of  the  analysis,  a rectangular  tank  is  assumed  of 
wetted  and  unwetted  wall  heights  H^  and  H^,  width  B and  length  W (as  shown 
in  Figure  II- 3).  The  effective  vent  area  is  A,  with  P-V  valve  open. 

For  a fully  loaded  tank,  la  less  than  4Z  of  the  total  tank  depth. 

Also,  unless  several  tanks  are  empty  and  the  vessel  displacement  is  unusually 
low,  H^  is  rarely  more  than  several  feet  above  waterline.  The  range  of 
parameters  used  in  the  present  study  is  representative  of  chemical  tanker 
operation: 

Parameter : 

Representative  Value: 

Range  of  Values: 

__ 

Approximately  equal  to  6"  pipe  diameter. 


H 

V 

4 ft 

2 - 31  ft 


«L 

3.3  ft 
3.3  - 16.5 


H 

25.5  ft 
25.5  ft 


2* 

.215  ft 

.054  - .54  ft2 
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For  the  baseline  vapor  space  height  of  4 ft,  the  total  voltnne  of  the  vapor 

space  assumed  was  2,700  ft  with  B assumed  to  be  26  ft. 

3.  Liquid  Cargo  Properties 

The  latent  heat  of  vaporization  of  the  cargo  controls  the  vapor  genera- 
tion rate,  once  the  wall  was  heated  up  and  the  heat  flux  to  the  liquid  has 
equilibrated.  The  values  assumed  for  this  key  parameter  and  for  the  vapor 
molecular  weight  are  given  below. 

Heat  of  Molecular 

Parameter : Vaporization  (X  ) Weight  (yv) 

Representative  Value:  145  Btu/lb  66  Ib/lb  mole 

Range  of  Values:  86  - 259  Btu/lb  20  - 100  Ib/lb  mole 

Representative  values  correspond  to  Hexane. 

4.  Failure  Pressure  of  Unwetted  Tank  Wall 

One  can  use  plastic  deformation  theory  to  estimate  the  failure  pressure 
at  room  temperature.  (See  Vent  System  and  Loading  Criteria  for  Avoiding  Tank 
Overpressurization,  part  3 of  Contract  DOT-CG-42,357-A,  p.  14  ff.)  Using  a 
modified  theory  of  Clarkson  (1956),  estimates  were  made  of  the  minimum  Internal 
pressure  loadings  required  to  Initiate  failure  of  cargo  tanks  of  three  representa- 
tive vessel  designs.  For  the  offshore  barge,  an  estimate  of  6 pslg  was  obtained; 
for  the  Inland  barge,  4 pslg;  and  for  a large  tankshlp,  8 pslg.  The  average  value 
for  these  three,  6 pslg,  will  be  utilized  In  this  report  as  a nominal  failure 
pressure  for  tank  vessel. 

As  the  wall  temperature  Increases,  the  material  strength  decreases.  The 
strength  does  not  decrease  substantially  below  about  600F,  as  shown  In  Figure  II-4. 
The  strength-temperature  relationship  appears  to  be  approximately  log  linear  between 
900F  and  1400F  for  steel,  with  a decrease  of  about  25Z  for  each  lOOF  above  900F.  The 
data  In  Figure  II-4 

* The  failure  pressure  for  a particular  tank  may  be  less  or  greater  than  this  value, 
and  will  depend  on  the  structure  and  condition  of  the  tank 
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for  two  structural  steels  is  taken  from  Weiss,  et  al.  (1966).  If  we  assume 
that  Clarkson's  equation  holds  and  P is  proportional  to  ^ /a  , then 

QAX 

Figure  II-4  can  be  used  to  estimate  the  decreasing  ability  of  the  wall  to 
withstand  pressure  In  the  vapor  space.  For  example.  If  an  unwetted  wall 
normally  capable  of  6 pslg  should  reach  1400F  temperature,  an  estimate  of 
the  failure  pressure  would  be  16Z  of  6 pslg  or  1.0  pslg.  Here  the  16% 
factor  was  obtained  from  Figure  11-4  at  1400F. 

An  additional  stress  factor,  which  was  not  dealt  with  In  the  present 
analysis.  Is  the  steep  thermal  gradient  along  the  exposed  wall  both  at  the 
water  line  and  at  the  liquid  cargo  line. 
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111.  ANALYSIS  AND  MODEL  FORMULATION 


A.  Physical  Mechanisms 

For  the  purposes  of  discussion  and  calculations,  we  consider  a 
rectangular  tank,  one  wall  of  which  is  exposed  to  the  fire.  The  physical 
system  is  illustrated  schematically  in  Figure  III-l,  and  the  Important 


FIGURE  MI-2  PHYSICAL  PROCESSES  DETERMINING  THE  PRESSURE  RISE 
DUE  TO  EXTERNAL  HEATING 


The  radiation  from  an  external  fire  heats  up  the  portion  of  the  wall  which 
is  above  the  water  line.  That  part  of  the  wall  in  contact  with  the  vapor 
heats  up  more  rapidly  than  the  one  in  contact  with  the  liquid.  The  wall  In 
contact  with  the  liquid  transfers  heat  to  the  liquid,  thereby  setting  a 
flow  of  the  heated  liquid  near  the  wall.  The  difference  between  the 
heat  coming  from  the  fire  into  the  wall  and  the  heat  expended  In  heating 
the  fluids  (vapor  or  liquid)  on  the  inner  face  goes  to  increase  the  wall 
temperature. 

When  a liquid  in  a tank  is  heated  from  the  sldea,  it  Is  well  known 
that  the  resulting  flow  of  liquid  produces  a density  stratification  pheno- 
menon irtilch  is  very  stable.  In  addition,  for  reasonably  high 
haat  flux  input  into  the  wall,  nuclaate  boiling  takes  place. 
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It  has  been  found  that  in  spite  of  the  liquid  boiling  on  the  hot  wall 
surface,  the  stratification  of  the  bulk  of  liquid  persists  [Anderson, 
et  al.,  (1974)].  Gebhart  (1971)  has  reported  that  up  to  10^  Btu/hr/ft^ 
heat  flux  can  be  accepted  from  a submerged  heater  by  nucleate  boiling. 

The  mass  of  vapor  emanating  from  the  liquid  will  Increase  substantially 
when  boiling  occurs  on  the  wall. 

Once  the  liquid  starts  boiling  on  the  wall  (and  even  if  the  bulk 
of  liquid  is  not  at  saturation  temperature) , the  temperature  of  the  wall 
in  contact  with  the  liquid  remains  essentially  unaltered  because  of  the 
very  high  heat  fluxes  that  occur  in  the  boiling  process.  The  average 
wall  temperature  will  be  between  10  and  30F  higher  than  the  saturation 

temperature  of  the  liquid.  Bolling  results  in  considerably  increased 
vapor  flow  into  the  ullage  volume.  Consequently,  the  tank  pressure  builds 
up  rapidly.  In  the  ullage  volume,  the  confined  vapor  is  heated  directly 
by  the  hot  unwetted  wall  by  convection.  Also,  mixing  of  the  vapor  generated 
by  boiling  with  the  pre-existing  ullage  contents  occurs,  and  molecular 
weight  differences  must  be  considered.  The  Increased  pressure,  together 
with  diminishing  strength  of  the  wall  in  contact  with  the  vapor  due  to 
Increased  temperature,  may  result  in  situations  where  the  tank  would  fall. 
The  schematic  diagram  III-2  illustrates  these  basic  physical  phenomena. 

In  this  report,  we  have  analyzed,  by  suitable  models,  the  above 
physical  phenomena.  The  pressure-time  history  and  the  wall  temperature- 
time  history  are  calculated  for  a variety  of  tank,  liquid,  fire,  and 
vent  characteristics. 

B.  Model  Formulation  and  Assumptions 

The  model  developed  to  determine  the  pressure-time  history  Involves 
essentially  writing  the  appropriate  mass,  energy,  and  specie  balance 
equations  in  addition  to  accounting  for  the  vapor  heating  and  liquid 
boiling.  Figure  II1-3  Illustrates  schematically  the  various  heat  sinks 
in  the  system  that  absorb  the  heat  from  fire.  We  write  the  equations 
to  describe  the  state  of  the  system  subject  to  the  following  assumptions: 
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Convective  Heating 
of  Vapor 


Heat  Radiated 
From  Inner  Wall 


Wall 

Temperature 

Constant 


Liquid  Boils 
on  the  Wall 


X Is  Wall  \ 

/ Temperature 
Greater  Than  Liquid 
Saturation  > 
\T  emperature7y^ 


Liquid  Stratifies 
Very  Little 
Evaporation 


Part  Reaches  Liquid, 
Resulting  in 
Evaporation  at  Surface 


FIGURE  III-3  HEAT  BALANCE  DIAGRAM  WHICH  SCHEMATICALLY  INDICATES 
THE  VARIOUS  SINKS  FOR  THE  HEAT  FROM  THE  FIRE 


Assumptions 


• The  liquid  is  characterized  by  a single  boiling  point,  and  its  | 

vapor  has  unique  thermodynamic  properties.  ! 

• The  vapor  produced  and  any  air  present  in  the  ullage  volume  act 
like  perfect  gases. 

• Mixing  of  vapor  and  air  in  the  ullage  occurs  Instantaneously  so 
that  there  are  no  temperature  or  composition  gradients  in  the 
mixture . 

i 

• The  change  in  the  liquid  volume  due  to  heating  is  small;  that  | 

is,  the  volume  of  the  ullage  space  remains  a constant.  No  1 

1 

cargo  loading  or  unloading  is  assumed  to  occur  during  the  fire.  | 

• Quasi-steady  analysis  is  used  to  describe  the  free  convection  ! 

(turbulent)  heat  transfer  between  the  wall  and  the  cargo  or  ! 

vapor;  that  is,  the  boundary  layer  is  assumed  to  be  fully 

developed  and  assumed  to  correspond  at  every  Instant  of  time 

to  the  wall  temperature  and  fluid  temperature  values  at  the  ! 

Instant  under  consideration. 

• The  radiation  from  the  wall  is  not  absorbed  by  the  gas  in  the 

ullage.  * 

1 

• The  heat  flux  from  the  fire  is  uniform  over  the  entire  wall  face  < 

i 

exposed  to  the  fire  and  constant  in  time.  | | 


• There  are  no  temperature  variations  across  the  thickness  of  i 

! 

the  wall;  that  is,  there  is  no  spatial  variation  of  temperature  i 

in  either  the  wetted  or  unwetted  wall.  However,  the  temperature 
is  a function  of  time. 

• The  temperature  of  the  wall  in  contact  with  liquid  remains 

constant  when  the  liquid  starts  boiling.  Nucleate  boiling  { 

rather  than  transition  or  film  boiling  is  assumed  to  occur.  | 

I 

• Heat  flux  reaching  the  cargo  surface  by  unwetted  wall  reradiation 
is  assumed  to  generate  vapor  rather  than  heat  up  the  liquid. 

a In  spite  of  the  large  temperature  gradient  in  the  wall  near 
the  liquid  meniscus,  the  total  heat  flow  rate  by  conduction 
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along  the  wall  Is  negligible.  Alternatively,  the  temperatures 
of  the  wair  . are  not  coupled  by  heat  transfer  between  them.  Con- 
duction to  the  cooler  deck  or  supporting  structure  Is  also 
neglected. 

• Radiative  heat  received  by  colder  walls  In  the  tank  Is  not 
reradlated  to  the  liquid. 

• The  cargo  Is  carried  at  ambient  pressure  at  the  start  of  the  fire. 

Hodel  Equation  Formulation 

We  write  below  the  various  "balance"  equations.  The  definitions  of 
the  symbols  are  given  In  the  nomenclature.  The  various  terms  are  also 
defined  schematically  in  Figure  III-4. 

C.  Governing  Equations 

1.  Thermal  Balance  on  the  Steel  Wall  In  Contact  with  Vapor 


dT 

P c L 
w w w dt 

I I 


Pate  of  Increase  in  the 
Internal  energy  of  a 
unit  area  of  wall 


Net  heat  coming 
into  the  wall 
per  unit  area 


where,  assuming  no  reradiation  into  the  wall. 


(1) 


(2) 


and  q"  is  the  mean  heat  flux  leaving  the  wall  due  to  turbulent  free  con- 
c ~ 

vection  cooling  by  the  vapor  flow.  The  expression  for  q'^  is  obtained  from 

the  expression  given  in  Appendix  B (equation  B-4) . 


0.0246 


K (T  - T)  rgB(T  - T)  H 

w w I 

H 2 

V L V 


, 2/5 

-] 


Pr 


7/15 


n 


2/3 

+ 0.494  Pr  '^1 


2/5 


(3) 


It  Is  noted  that  In  the  above  expressions,  both  T and  T (the  vapor  tem- 

w 

perature)  are  functions  of  time. 
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FIGURE  III-4  ILLUSTRATION  OF  HEAT  TRANSFER  TERMS 


2.  Heating  of  Liquid  and  Its  Subsequent  Boiling 

When  the  wall  wetted  by  the  liquid  Is  exposed  to  radiation  from 
fire,  a part  of  the  heat  is  absorbed  by  the  wall  material  and  the  re- 
maining part  goes  to  heat  the  liquid  layers  adjacent  to  the  wall. 

Heating  of  the  liquid  creates  a free  convection  flow.  The  heated  liquid 
accumulates  in  the  top  layers  of  the  liquid  bulk  in  the  form  of  a stably 
stratified  layer.  The  depth  of  this  layer  increases  continuously  with 
time.  An  analysis  of  heating  of  the  liquid  which  considers  this  liquid 
stratification  phenomenon  is  developed  in  Appendix  A.  It  is  found  from 
this  analysis  that  for  a typical  tank  and  liquid  combination  representing 
the  most  common  types  of  tanks  and  liquids  used  in  practice,  the  time  of 
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exposure  to  typical  flame  radiation  (32,000  Btu/hr  ft  ) necessary  before 
the  wall  temperature  attains  boiling  temperature  of  liquid  Is  very  small 
(of  the  order  of  100  seconds). 

Once  the  wall  temperature  reaches  a few  degrees  above  the  liquid 
boiling  temperature.  It  Is  reasonable  to  expect  the  liquid  to  boll.  Since 
the  time  for  the  wall  to  attain  the  boiling  temperature  Is  small,  we 
assume,  for  pressure  calculation  purposes,  that  boiling  ensues  from  the 
Instant  the  wall  Is  exposed  to  fire  radiation.  The  result  of  this  assump- 
tion would  be  to  predict,  during  the  initial  period,  a quicker  rate  of 
Increase  of  pressure  than  that  which  would  occur  If  the  transient  heating 
of  and  stratification  of  liquid  phenomena  are  taken  Into  account.  There- 
fore, using  the  equations  given  below  would  amount  to  a conservative 
prediction  of  the  time  of  tank  rupture. 

Since  the  heat  flux  to  the  liquid  during  boiling  Is  very  much 
larger  than  when  the  liquid  is  being  heated  convectively , the  wall 
temperature  remains  a constant  once  boiling  ensues.  Therefore, 


(4a) 

(4b) 


We  further  assume  that  all  of  this  vapor  generated  by  the  liquid  boiling 
Is  released  Into  the  vapor  space  In  the  tank;  that  Is,  no  part  of  vapor 
dissolves  during  Its  ascent  through  the  liquid. 


Additional  mass  of  vapor  Is  released  Into  the  ullage  volume  because 
of  the  evaporation  of  liquid  at  Its  surface.  The  heat  for  this  evapora- 
tion comes  from  the  hot  wall  In  contact  with  the  vapor  by  direct  radiation. 


M - W H 
evap  v 


F q" 
rad 


(5) 


where  F la  the  view  factor  for  radiative  exchange  between  the  hot  wall 
in  contact  with  the  vapor  and  the  liquid  surface  (which  la  assumed  to 
be  flat).  The  expression  for  calculating  F has  been  taken  from  Hiebelt  (1966) 
and  is  given  In  Appendix  B.  Therefore, 


■ M.  + M • total  rate  at  which  vanor  Is 
D evap 


■’Injected’’ 
Into  the  ullage  volume 
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Global  Mass  Balance  Equation 


We  now  write  the  mass  balance  equations  for  air  and  vapor  species 
as  well  as  the  total  mass  of  gases  within  the  ullage  volume. 
dM 

dt  “ “l  - ”o  (7) 

• 1 1 I 

Rate  of  change  of  mass  Net  Inflow  of  gases 

In  the  ullage  volume 

Specie  Mass  Balance  Equations 
d(Y  M) 

■ ■ ^a  ^ (8a) 

Rate  of  change  of  Mass  outflow  rate 

mass  of  air  In  the  of  air 

ullage  volume 


Similarly,  the  mass  balance  equation  on  the  vapor  Is  written  as 


d(Y^M) 

dt~ 


Y M 


and  It  Is  noted  that 


(Sb) 


1 « Y + Y 
a v 


(9) 


In  writing  equations  8a  and  8b,  It  Is  assumed  that  no  air  Is  generated 
by  the  boiling  process  and  that  In  the  venting  gas,  the  proportion  of 
air  and  vapor  (by  mass)  is  the  same  as  In  the  ullage  volume. 

Using  equations  7 and  8a,  we  can  write  the  air  mass  balance 
equation  as 


Venting  Equation 

According  to  Bernoulli's  law,  the  mass  flux  through  a duct  of  cross 
sectional  area  A can  be  related  to  the  pressure  drop: 

K ■ * ' 

where  p is  the  density  of  gases  In  the  ullage  space  and  C.  is  the  dls- 

it  ^ 

charge  coefficient  . 


Note  that  the  pressure  difference  Is  In  psi,  the  area  In  ft^,  and  the 
density  Ibm/ft^  so  that  a conversion  factor  of  12  In/ft^  must  be 
applied  to  the  right  hand  side  when  using  equation  (11). 
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Pressure  Equation 

The  pressure  In  the  gas  space  is  due  to  the  partial  pressures  of 
air  and  vapor. 


with  p 


p - + 

Pv 

(12a) 

M 

R T 

a 

P«  “ TT* 

u 

(12b) 

‘^a 

M 

R T 

V 

Pv“  v" 

u 

w 

V 

(12c) 

A, 

dM  T 

dM  > 

k /M  M > 

dT 

— - + — 
dt 

V 

+ ( — + — 

(12d) 

Differentiating  12a  with  respect  to  time,  and  using  12b  and  12c,  we  get 


dt 


Energy  Conservation  Equation 

In  order  to  estimate  the  temperature  variation  of  the  gas  with 
respect  to  time  when  It  is  subjected  to  heat  from  the  walls  and  when 
there  is  venting,  we  write  the  energy  equation.  The  derivation  of 
the  equation  Is  Involved  and  quite  extensive;  therefore.  It  Is 
given  in  Appendix  C.  Only  the  final  result  Is  Indicated  below. 


dt 


"r  *1 

r Y Y 1 

• • 

— T + cP(T.  - T) 

- M 

— R T 

Q M 

y v'  b 

o 

y y„  u 

c 1 

, 

.a  V . 

V V 

M[Y  c + Y c'^l 
a a V V 


(13) 


Equations  (1),  (10),  (12d) , and  (13)  represent  the  variation  with  time  of 
the  four  unknowns  in  the  system;  namely,  T^,  Y^,  p,  and  T.  All  other  equa- 
tions contain  parameters  that  can  be  expressed  in  terms  of  the  above  unknowns. 
The  problem  therefore  reduces  to  solving  for  four  unknowns  from  four  simul- 
taneous, coupled,  first  order,  nonlinear  differential  equations.  The  solu- 
tion procedure  is  briefly  explained  In  the  next  section.  The  initial  values 
of  the  paraawters  used  depend  on  the  type  and  dimensions  of  the  tank, 
amount  of  air  In  the  ullage  volume,  etc. 
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D.  Solution  Procedure 


The  solution  to  the  above  equations  Is  obtained  by  first  writing 
all  of  the  equations  In  dimensionless  form.  This  Is  Illustrated  In 
Appendix  D.  Then  the  nondlmenslonallzed  equations  1,  10,  12d,  and  13 
are  solved  numerically  on  a computer  using  the  Hamming  Predictor  Corrector 
Technique  [Hamming  (1973)].  This  technique  Is  essentially  a procedure 
which  Involves  Incrementing  time  In  a stepwise  fashion  finding  the 
rate  of  change  of  parameters  at  the  given  time  step  (with  known  values 
at  the  given  time).  Then  the  values  are  calculated  at  the  next  time 
step  and  errors  evaluated  between  the  calculated  values  and  the  values 
obtained  If  one  half  of  the  time  step  Is  used.  If  the  error  Is  greater 
than  a specified  value,  the  time  step  Is  halved  and  the  procedure 
repeated.  The  program  Is  capable  of  halving  the  time  step  ten  times 
automatically  to  achieve  the  desired  error  bound  (l.e.,  the  lowest 
values  of  time  step  that  the  program  will  use  will  be  equal  to  the 
given  Initial  time  step  divided  by  2^®) . 


E.  STEADY  STATE  SOLUTION 


If  the  fire  continues  Indefinitely,  the  tank  pressure  vapor  temperature, 
and  unwetted  wall  temperature  will  gradually  level  off  at  steady  state  val- 
ues. Under  these  conditions,  the  air  Is  completely  displaced  by  vapor  so 
that  p » u^,  - 0,  and  - 1.  In  addition,  the  unwetted  wall  tempera- 

ture Is  quite  large  (above  800C  according  to  results  to  be  presented  In 
Chapter  IV),  so  that  convection  can  be  neglected  relative  to  reradlatlon 

(q^"  ■ 0)  . Applying  these  assumptions  and  setting  time  derivatives  equal 
to  zero  In  Equations  (7)  and  (13),  one  can  obtain  the  following  solution  for 
T and  p: 

T - T. 

D 

Using  Equation  (12c)  to  express  p In  terms  of  p and  T,  and  substituting 
(14)  Into  (15),  we  can  derive  a quadratic  equation  for  pressure: 


(14) 

(15) 


The  solution  to  Equation  (17)  Is  as  follows: 


^ **amb  ^arab 


y/1  + - l] 

L 2 J 


Calculations  show  that  at  a wall  temperature  of  730C,  q " ■ 6Z  of  q,",  com- 
pared to  4 ■ 94Z  of  4.”,  if  the  emlssivlty  Is  unlty^(as  assumed‘ln 

Equation  (2)).  * 


where  C 


2“vVp..b‘*'=dV  ^// 


Under  the  condition  C l/4>  Equation  (18)  can  be  simplified  to  read 


P - P V " P w 

amb  amb 


(Large  C)  (19) 


However,  for  baseline  values  (see  Chapter  IV),  C ==  .05  and  so  the  assump- 
tion 4 C <<  1 more  applicable.  Under  this  assumption. 


P-P  u=P  b.5- 

amb  amb 


(Small  O (20) 


In  the  regime  of  parameter  values  for  which  Equation  (20)  Is  valid,  the 
steady  state  pressure  in  the  tank  is  predicted  to  vary  as  follows: 


P-  P , 
amb 


R T,. 
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2m  p ,\AAC. 
V amb  ' d 
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(21) 
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IV.  RESULTS  AND  DISCUSSION 


A.  PRESSLRL  RISE  AND  UNVETTED-WALL  TEMPERATliRL  BASELINE  CASE 


The  computer  simulation  described  in  Part  III  for  the  cargo  tank  ex- 
posed to  external  fire  was  exercised  for  the  following  baseline  values: 


Fire  heat  flux 
Cargo  latent  heat 
Cargo  molecular  weight 
Wetted  wall  area 
Unwetted  wall  height 
Effective  vent  area 


(32,000  Btu/hr/ft^) 

145  Btu/lb  ) (hexane) 

86  Ib/lb  mole 
84  ft^ 

4.  ft 

.215  ft^  (6"  vent  valve) 


The  pressure  rise  under  these  conditions  levels  off  at  about  1.5  psig  in 
about  25  minutes,  as  shown  in  Figure  IV-1.  The  relatively  low  level  of 
pressure  is  significant;  to  the  extent  that  the  input  data  is  realistic 
this  means  that  many  cargo  tanks  are  readily  capable  of  venting  the  vapor 
generated  by  fire.  The  pressure  rise  for  a lighter  vapor  fraction  (u^  ■ 30, 
say)  could  be  as  high  as  3.5  psig  as  will  be  shown  below. 

An  Interesting  qualitative  feature  of  the  pressure  rise  is  that  it 
occurs  in  two  stages.  A rapid  Increase  to  about  0.5  psig  within  20 
seconds  occurs  due  to  Immediate  boiling  from  the  wetted  wall.  A more 
gradual  pressure  rise  follows,  caused  by  the  heat  up  of  the  unwetted 
wall  (which  in  turn  heats  up  the  confined  vapor).  The  indentation  notl- 
cable  in  Figure  IV-1  at  about  500  sec  is  pressumably  ascrlbable  to  the 
gradual  shift  in  average  molecular  weight  from  the  original  tank  contents 
(assumed  * 30)  to  the  pure  vapor  (w^  ■ 86). 

That  the  venting  system  can  accommodate  the  vapor  generated  by  a 
2 

32,000  Btu/hr/ft  external  fire  with  only  a 1.5  psig  pressure  rise  does  not 
guarantee  that  the  tank  will  remain  intact,  of  course.  Figure  IV-2 
shows  that  the  unwetted  wall  temperature  exceeds  600C  in  about  10  minutes. 

In  15  minutes,  the  wall  strength  will  have  been  reduced  by  heating  to  the 
point  where  failure  at  2 psig  is  entirely  possible  for  a representative 
tank  (normally  capable  of  6 psig).  This  result  suggests  that  the  failure 
node  of  cargo  tanks  is  due  not  to  the  vent  system,  but  due  to  simple 
weakening  of  the  umratted  wall.  It  also  suggests  that  the  method  of 
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analysis  of  unwetted  wall  temperature  should  be  refined  by  Including 

lateral  heat  conduction  to  supporting  structure,  stiffeners,  and  the  deck 

i plating  not  exposed  to  fire.  The  key  Input  data  In  our  prediction  of 

800C  wall  temperature  can  also  be  scrutinized:  Fire  flux  of  32,000 
2 

Btu/hr/ft  , unwetted  wall  height  of  4 ft,  and  plating  thickness  of 


.062  ft  (.75  In). 


B.  EFFECT  OF  HEAT  FLUX  LEVEL 


2 

The  heat  flux  level  reaches  32,000  Btu/hr/ft  only  in  the  extreme  case  of 

a large  fire  Immediately  adjacent  to  the  cargo  tank.  If  the  heat  flux  level 

is  less  than  this,  the  unwetted  wall  temperature  will  be  less  and  also  the 

vapor  generation  rate  will  be  less.  In  Figure  IV-3  we  have  ;>hown  predictions 

of  the  wall  temperature  for  heat  flux  ranging  from  16,000  to  64,000  Btu/hr/ft  . 

2 

Note  that  at  16,000  Btu/hr/ft  the  wall  temperature  does  not  exceed  S50C  in 

the  first  20  minutes  of  the  fire,  and  making  reference  to  Figure  II-4  one  can 

deduce  that  the  tank  pressure  can  rise  to  60%  of  the  normal  maximum  without 

2 

failure.  On  the  other  hand  extreme  fire  fluxes  of  150,000  Btu/hr/ft  are 
predicted  to  cause  local  wall  collapse  within  10  minutes  (under  the  assump- 
tions of  the  present  model). 

The  pressure  rise  at  elevated  heat  flux  is  displayed  in  Figure  IV-4. 

2 

If  the  heat  flux  doubles  from  32,000  to >64,000  Btu/hr/ft  , the  pressure  rise 

after  20  minutes  apparently  increases  about  a factor  of  four.  The  re- 

• 2 

latlonshlp  AP  q can  be  explained  on  the  basis  that  AP  is  proportional 
to  the  square  of  the  velocity  in  the  vent,  which  is  in  turn  proportional 
to  the  vapor  generation  rate  caused  by  the  heat  flux.  For  a 6"  vent 
diameter,  the  pressure  rise  is  less  than  6 pslg  even  with  the  maximum 
conceivable  heat  flux  (64,000  Btu/hr/ft  ).  This  reiterates  the  earlier  con- 
clusion that  the  external  fire  hazard  is  attributable  to  the  unwetted 
wall  rather  than  the  vent  system. 
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C.  EFFECT  OF  CARGO  HEAT  OF  VAPORIZATION 


In  steady  state  after  wall  heat  up,  the  rate  of  vapor  generation  Is 
Inversely  proportional  to  the  heat  absorbed  per  unit  vapor,  the  largest 
part  of  which  is  the  heat  of  vaporization.  Thus,  the  velocity  of  vent 
gases  should  be  nearly  Inversely  proportional  to  the  heat  of  vaporization. 
Since  AP  v^,  one  therefore  expects  AP  'v  A 

Since  most  hydrocarbons  show  a heat  of  vaporization  in  the  range  100 
to  300  Btu/lb...  we  have  selected  hexane  (145  Btu/lb)  for  the  baseline 
representative  value.  We  have  compared  the  pressure  rises  for  several 
values  of  vaporization  heat  in  Figure  IV-5.  For  the  lowest  reasonable 
value  (86  Btu/lb),  the  pressure  rise  is  still  less  than  3.2  psig  after 
20  minutes,  which  is  quite  safe  for  a wall  of  normal  strength. 
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FIGURE  IV-6  EFFECT  OF  HEAT  OF  VAPORIZATION 
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D.  EFFECT  OF  EXPOSED  CARGO-TANK-WALL  AREA 


In  the  usual  case,  the  amount  of  liquid  cargo  stored  above  the  water- 
line is  not  large;  in  fact  we  have  assumed  an  exposed  wetted  wall  area 
2 

of  only  84  ft  per  tank  wall.  If  several  adjacent  tanks  are  empty  of  for 
some  other  reason  the  tank  exposed  to  fire  protrudes  well  above  the  water- 
line, more  vapor  will  be  generated  by  the  fire.  Figure  IV-6  shows  how 
the  exposed  wetted-wall  area  affects  the  pressure  rise.  A factor  of  three 
Increase  in  area  leads  to  a factor  of  four  Increase  in  pressure  differential 

after  20  minutes  of  fire  exposure.  If  the  exposed  wetted  area  climbs  to 
2 

420  ft  (16  feet  above  waterline  for  a 26  ft  wide  tank),  the  pressure  is 
predicted  to  rise  to  values  exceeding  6 pslg.  In  view  of  the  heat  up 
and  weakening  of  the  unwetted  wall,  this  would  be  an  extremely  hazardous 
situation. 

2 

Suppose  on  top  of  the  abnormally  high  (420  ft  ) exposed  area,  one 
assumes  that  the  cargo  liquid  has  a reasonably  low  heat  of  vaporization 
(86  Btu/lb).  Then,  as  shown  in  Figure  IV-7,  the  pressure  rise  is  predicted 
to  exceed  6 pslg  and  is  almost  certain  to  cause  failure  even  for  unwetted 
walls  of  normal  strength. 
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FIGURE  IV-7  EFFECT  OF  COMBINED  SMALL  HEAT  OF  VAPORIZATION 
AND  LARGE  EXPOSED  AREA 
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E.  EFFECT  OF  ULLAGE  VOLUME 


The  ullage,  or  portion  of  the  tank  not  filled  with  liquid  cargo,  in- 
creases the  pressure  rise  by  transmitting  additional  heat  to  the  confined 
vapor  from  the  inside  of  the  unvetted  wall.  This  occurs  by  convection 

to  the  vapor  and  by  radiation  to  the  liquid  surface.  A series  of  runs  was 

2 

conducted  at  32,000  Btu/hr/ft  heat  flux,  assuming  Hexane,  varying  only 
the  height  of  the  vapor  space  from  2 to  32  feet  (baseline  4 feet). 

The  ullage  volume  only  affects  the  later  stages  of  pressure  rise,  as 
shown  in  Figure  IV-8.  At  7 minutes  from  the  start  of  fire  exposure,  the 
pressure  rise  is  0.7  psig  for  all  cases  and  is  mainly  due  to  vapor  gen- 
eration (not  "afterheating"  by  the  unwetted  wall).  At  later  stages, 
say  after  20  minutes,  the  unwetted  wall  is  transferring  more  heat  to  the 
confined  vapor  and  the  pressure  now  is  quite  sensitive  to  ullage  volume. 

A factor  of  8 increase  in  vapor  space  leads  to  an  Increase  in  tank  pres- 
sure from  1.3  psig  to  3.0  psig  (both  after  20  minutes). 

The  sensitivity  of  pressure  to  ullage  volume  is  not  so  great  that 
the  venting  capacity  becomes  limiting;  the  maximum  pressure  in  the  cases 
considered  was  3.0  psig. 
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FIGURE  IV-8  EFFECT  OF  PARTIALLY  EMPTY  TANK 


F.  EFFECT  OF  VENT  AREA 


Vent  pipes  are  normally  sized  according  to  the  cargo  Inlet  pipes, 

which  In  turn  are  sized  according  to  the  tank  valves  so  that  the  tank  can 

be  filled  In  4 - 8 hours  with  about  10  ft/sec  liquid  velocity.  The  value 
2 

of  .215  ft  adopted  for  the  baseline  vent  area  should  correspond  to  the 
inlet  pipe  size  on  a 82  x 26  x 26  ft  tank  since  this  gives  25,800  seconds 
loading  time  (7.2  hours),  at  10  ft/sec  loading  velocity. 

Smaller  tanks,  say  on  barges,  would  have  smaller  vent  areas,  down  to 
2 

.054  ft  (approximately  3"  pipe),  and  still  might  have  the  exposed  wetted 
2 

area  of  84  ft  we  are  using  as  baseline.  The  pressure  rise  is  obviously 
much  larger  for  such  cases,  as  shown  In  Figure  IV-9.  The  pressure  rise 
may  be  as  great  as  10  pslg  for  a .054  £t^  vent  area  according  to  our  pre- 
dictions. In  practice,  a 3"  vent  would  most  often  be  found  on  a smaller 

2 

tank  with  much  less  than  84  ft  exposed  wetted  wall,  and  the  pressure  would 
be  correspondingly  lower.  In  general,  our  results  suggest  that  the  exposed 
wall  area  should  not  be  more  than  1,000  times  the  vent  area  If  the  pressure 
Is  not  to  exceed  6 pslg  for  a typical  cargo  and  heat  flux. 


Pressure  (psig) 


FIGURE  IV-9  EFFECT  OF  VENT  AREA 
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G.  EFFECT  OF  VAPOR  MOLECULAR  WEIGHT 


Since  the  volume  taken  up  by  a gas  is  inversely  proportional  to  its 
molecular  weight,  the  confined  vapor  generated  by  an  external  fire 
causes  a greater  pressure  rise  the  lower  the  molecular  weight.  This  is 
Illustrated  by  Figure  IV-10;  the  pressure  rise  after  20  minutes  is  3.4 
pslg  for  a 30  Ib/lb  mole  vapor  (compared  to  a baseline  of  1.3  pslg  for 
hexane  at  86  Ib/lb  mole) . The  range  of  molecular  weights  has  a lower 
limit  of  16  Ib/lb  mole  (CH^)  so  that  hazardous  pressures  are  not  expected 
due  to  the  molecular  weight  variation. 

No  physical  significance  should  be  attached  to  the  "spikes"  which 
appear  in  Figure  IV-10  at  vapor  molecular  weights  of  50  and  30.  Physi- 
cally, at  the  onset  of  a fire  there  is  almost  Immediate  boiling  at  the 
wetted  wall.  Numerically,  the  venting  rate  is  initially  zero  and  under- 
predicted  for  the  first  few  time  steps.  This  results  in  an  over- 
prediction  of  dP/dt  (see  Equation  (12d)]  initially  which  is  magnified 
for  low-molecular-weight  vapor  because  the  is  in  the  demoninator. 

The  magnitude  of  the  steady  state  pressure  rise  can  be  shown  to  be 
inversely  proportional  to  the  vapor  molecular  weight  for  small  pressure 
rises,  in  agreement  with  Figure  IV-10.  Equation  (E-6)  would  take  the 
following  form  for  variable  mixture  molecular  weight: 


yi  dP 

Wa  dt 


a 


Under  steady  state  conditions  the  left  hand  side  is  zero,  p > and 

the  dimensionless  state  pressure  for  small  pressure  size  is  given  by 
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H.  COMBINED  EFFECT  OF  SEVERAL  PARAMETERS  AND  HAZARD  EVALUATION  CRITERION 
1.  Correlation  of  Predicted  Data 

The  results  presented  In  Sections  A through  G describe  the  effect  of 

each  single  variable  on  the  pressure  rise.  In  most  cases  the  pressure 

levels  off  after  about  20-30  minutes  to  a steady  state  value.  In  this 

section  we  relate  this  maximum  pressure  rise  p , to  arbitrary  combination 

ss 

of  the  parameters  q^,  X,  W,  and  A.  A closed  form  solution  of 

Equations  1,  10,  12d,  and  13  can  be  used  to  derive  a general  pressure  rise 
equation  applicable  to  any  combination  of  variables. 


A hazard  criterion  covering  arbitrary  combinations  of  the  variables 


q^",  X,  W,  and  A can  be  obtained  from  Equation  (18),  which 

Is  the  asymptotic  value  of  tank  pressure  reached  after  about  30  minutes 
of  fire  exposure: 


amb 


^ss  ^amb  2 


where  we  have  taken  C , 


/ ^f"  ^ 
2 \j  ^ Pamb  ^ 

1 for  simplicity. 


If  the  tank  pressure  p - P , must  not  exceed  a certain  value,  then 

*^ss  amb 

Equation  (18)  can  be  shown  to  require  that 


(23) 


< 1.3  X 10 


10 


2 

where  q " Is  In  Btu/hr/ft  , X Is  In  Btu/lb,  and  p - p . Is  In  psl. 
r 88  amo 

Substituting  p^^  - P^jj  ■ 4 pslg,  we  obtain 


(22) 
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Equation  (24)  may  be  termed  the  4 pslg  criterion  for  venting  tanka 
exposed  to  external  fires.  It  should  be  noted  that  Equation  (22)  was 
derived  on  the  basis  that  convective  heat  transfer  was  negligible,  whereas 
the  numerical  results  of  Figures  IV-1  through  IV-10  Include  the  effect  of 
convection.  Therefore,  Equation  (22)  slightly  underestimates  the  pressure 
rise.  Sample  calculations  show  <1^,"/  to  be  from  5X  to  10%,  depending 
on  unwetted  wall  temperature,  emlssivlty  and  vapor  temperature.  Therefore, 
the  correction  to  maximum  pressure  due  to  convection  Is  not  expected  to 
be  significant. 


APPENDIX  A 


Thermal  Stratification  of  Liquid  Cargo 

Objective 

The  objective  of  this  appendix  is  to  formulate  a model  for  the 
stratification  phenomenon  that  occurs  when  a liquid  is  heated  and  to 
obtain  a time  history  for  the  temperature  of  the  wall  in  contact  with 
the  liquid. 

Physical  Model  j 

I 

The  physical  model  considered  is  shown  in  Figure  A-1.  The  heat  } 

from  the  fire  first  heats  the  wall  (in  contact  with  the  liquid)  which  | 

in  turn  heats  the  liquid  adjacent  to  it.  Because  the  heated  liquid 
becomes  buoyant,  it  rises  to  the  surface,  setting  up  convection. 

The  heated  liquid  rises  to  the  surface  of  the  liquid  bulk.  Because 

the  wall  temperature  Increases  continuously  with  time,  the  temperature 

of  the  liquid  in  the  convective  flow  adjacent  to  tlie  wall  also  Increases, 

and  this  liquid  accumulates  on  the  surface  of  the  liquid  bulk.  The 

successive  accumulation  of  hot  liquid  results  in  the  formation  of  a 

stable  stratified  layer  as  shown  in  Figure  A-1.  The  hottest  liquid 

is  at  the  top.  The  bulk  of  the  liquid  is  not  affected  by  this  process 

for  a considerable  time.  The  depth  of  the  stratified  region  grows  with  : 

time. 

Asssumptlons  j 

I 

I 

The  above  process  is  modeled  below  with  the  following  assumptions:  | 

I 

• Turbulent  convection  in  the  fluid  is  quasi-steady;  that  is,  at 
every  Instant  the  convection  boundary  layer  is  fully  developed 
and  corresponds  to  the  wall  condition  existing  at  that  time. 

• Temperature  of  the  wall  in  contact  with  the  liquid  is  uniform  ’ 

in  spite  of  the  existence  of  the  stratified  layer;  that  is,  the  ! 

wall  temperature  predicted  is  a vertically  averaged  value. 

• Heat  flux  from  the  wall  into  the  stratified  liquid  layer  is 
equal  to  the  average  heat  flux  (in  the  turbulent  natural  con- 
vection cooled  region)  over  the  remaining  portion  of  the  wall 
wetted  by  the  fluid  above  the  water  line. 
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FIGURE  A-1  ILLUSTRATION  OF  FORMATION  OF  STRATIFIED  LAYER 


Water  Line 


• The  turbulent  boundary  layer  begins  at  the  water  line. 

• The  bulk  temperature  in  the  nonstratlf led  liquid  does  not  change. 


I 

[ 


• Evaporation  (and  therefore  mass  loss)  from  the  liquid  surface 
Is  negligible. 

• Density  change  of  liquid  due  to  heating  is  small  and  Is  con- 
sidered Insignificant  for  mass  balance  purposes. 

Governing  Equations 

1.  Mass  Balance  on  the  Stratified  Layer 

"L'h,  f ■ “o  <*-» 

where  A,  Is  the  surface  area  of  the  liquid  in  the  tank  and  M is  the 

L.  c 

mass  of  heated  liquid  flowing  Into  the  stratified  layer  by  convection. 

is  evaluated  from  the  knowledge  of  the  turbulent  boundary  layer 
characteristics.  This  Is  shown  in  Appendix  B.  Only  the  final  answer 
is  presented  here. 


g 6 (Hj^-S)  (T 


M * 0 . 1 W p,  V, 
c L L 


0.4 


Pr®^^^[l  + 0.494  Pr^^^ 


0.4 


(A-2) 


2.  Energy  Balance  on  the  Stratified  Layer 

In  writing  the  energy  equation,  we  assume  that  the  average  heat 
flux  value  (over  the  height  - S)  can  be  used  to  account  for  the  heat 
Into  liquid  over  the  height  H^.  It  is  noted  that  this  assumption  is 
If  S that  is,  when  the  stratified  layer  depth  has  grown 

to  an  extent  comparable  to  the  depth  of  liquid  exposed  to  fire.  When 
this  happens,  the  average  heat  flux  is  dominated  by  high  heat  flux 
values  even  though  in  reality  such  high  heat  fluxes  cannot  exist  due 
to  the  Influence  of  the  heated  layer  in  contact  with  the  wall.  How- 
ever, within  the  above  limitation  (that  S « Hj^),  we  write  the  following 
energy  flow  equation: 


A-3 


(A-3) 


S 

«- V'l 

*-  z-O 

where  T Is  the  temperature  within  the  stratified  layer. 

Inherent  in  equation  A-3  is  the  assumption  that  the  heat  capacity 
of  the  small  amount  of  liquid  In  the  turbulent  convective  boundary 
layer  is  negligible.  We  further  assume  that  the  temperature  distribution 
in  the  stratified  layer  is  self-similar.  The  accuracy  of  this  assump- 
tion has  been  verified  by  Clark  (1968).  Therefore,  defining 

1 

d(z/S)  (A-A) 

B / 

z/S-0 

from  experimental  data,  it  has  been  observed  by  Clark  (1968). 
that  the  above  energy  integral  is  essentially  constant  with  a value 
close  to  0.43.  Hence,  substituting  equation  A-4  in  A-3,  we  have 

VlVsf  - VI  - 


I « energy  integral 


I (. 


3.  Wall  Heat  Up  Equation 


dT 


L P c 

WWW 


w 


dt 


(A-6) 


where  q"  , Is  the  mean  heat  transfer  rate  by  natural  convection  into 
c ,L 

the  liquid.  Expressions  for  this  are  derived  in  Appendix  B in  terms 
of  the  wall  and  liquid  temperatures. 


Equations  A-1 , A-5,  and  A-6  form  a coupled  set  of  nonlinear  dif- 
ferential equations  in  three  unknowns;  namely,  S,  T^.  and  T^.  these 

can  be  determined  numerically  as  functions  of  tire  after  expressing 
• • 

M and  q"  , in  terms  of  T - T„  and  H,  - S as  indicated  in  the  next  section. 
C W S L 
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Solution  to  Equations 

Before  solving  the  above  set  of  equations.  It  Is  necessary  to  express 
them  in  dimensionless  parameters.  We  therefore  define  the  following: 


^ - — 

«L 


dimensionless  depth  of  stratified  layer 


PA  To  taken  to  raise  the  temperature  of 


ch 


XL  B 


the  liquid  In  the  tank  above  the  water 
line  by  a temperature  equal  to  the  initial 
bulk  temperature 


t/t  , 
ch 

(T/Tg  - 1) 


dimensionless  mean  heat  flux  into  liquid 
c'-'f  by  turbulent  convection 


q./q" 


Wp  c L H ratio  of  thermal  capacity  of  steel 
“to  that  of  the  liquid 
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M t , mass  of  liquid  that  can  accumulate  over  t 
B ch  ch  ^ 


L«L 


“ when  generated  at  a rate  o' 


ch 


mass  of  liquid  in  the  tank  above  the  water  line 
Using  the  above  definitions  and  equation  A-2,  we  write  equation  A-1  as 

e»-‘ 

di  w 

Equation  A-5  becomes 

' ST  <«  V ■ » 


and  equation  A-6  becomes 
d e 


di 


- 1 - q 


with  initial  conditions 
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(A-7a) 
(A- 7b) 

(A-7c) 

(A-7d) 

(.^-7e) 

(A-7f) 

(A-7g) 

(A-7h) 
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(A-11) 


Also,  it  can  be  shown  (using  the  expressions  for  derived  In  Appendix 
B)  that 

q - r 0 (1  - 
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where  r ■ 0.0246 


- 0.4 
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Substituting  for  q from  equation  A-12  In  equation  A-10,  we  have 

^ (0)  e (1  - 

^ (.  e^)  . 1 - r — ^1-:^ 


and  rewriting  equation  A-8,  we  have 


Jjr  [(1  - 


■[W  <" 


Also  adding  equations  A-9  and  A-10  and  Integrating  and  using  Initial 
conditions  A-11,  we  get 

(IC6  +u)6)*t 
s w 

Equations  A-14  and  A-15  form  a coupled  set  of  equations  In  unknowns 
u 6^  and  (1  - c)>  These  can  be  solved  numerically,  and  from  equation 
A-16,  the  dimensionless  liquid  surface  temperature  can  be  obtained. 

Specific  Example 

A specific  example  Is  worked  out  with  the  following  numerical 
values  to  Illustrate  the  procedure  of  calculation  as  well  as  to  obtain 
a feel  for  the  value  of  the  parameters  obtained. 


Tank  Dimensions 


Width 


W ■ 40  ft 


Height  of  liquid  level  u ■ jg  fj. 
above  the  waterline  l 


Breadth  of  tank 


B - 17  ft 


(A-12) 


(A-13) 


(A-14) 


(A-15) 


(A-16) 


A-6 


Density 

Specific  heat 

Kinematic  viscosity 

Prandtl  number 

Coefficient  of 
volume  expansion 


p_  - 56.2  Ibm/ft^ 

L 

Cj^  ■ 0.5  Btu/lbm  F 
V.  - 0.03875  ft^/hr 

PtL-  10 

6 » 4 X 10"^/“F 


Bulk  liquid  - 62*F 

temperature 

Steel  Wall  Properties 

3 

Density  * 500  Ibm/ft 

Specific  heat  c^  " 0.13  Btu/lbm  °F 

Thickness  L * 0.0625  ft  (3/4") 
w 

Radiation 

* * 2 

Heat  flux  from  fire  = 35,000  Btu/hr  ft 

Value  of  energy  Integral  1 >^  0.43 

' (see  equation  A-4) 

The  values  of  the  various  parameters  defined  in  equation  A-7  are: 

t . ■ 25640  seconds 
ch 


A - 10.378 

r - 2.554 

(1)  ■ 8.5  X 10 


Table  A-1  shows  the  results  obtained  by  solving  equations  A-14,  A-15, 
and  A-16.  The  parameters  have  been  expressed  in  dimensional  quantities. 

It  is  seen  that  for  the  example  considered,  the  wall  temperature  in 
contact  with  the  liquid  Increases  to  more  than  llO^C  in  about  100  seconds. 
Tills  Illustrates  the  rapid  temperature  rise  of  the  wall.  Also  by  noting 
that  the  stratification  depth  is  only  a few  Inches,  we  find  that  our 
assumption  S <<  H is  Justified.  Also  noteworthy  is  the  small  surface 

La 

temperature  rise.  This  further  Justifies  the  neglect  of  liquid  mass  loss 
by  surface  evaporation. 


* Calculation  was  completed  before  the  representative  value  of  32,000  Btu/ 
hr/ft^  (given  on  page  9)  was  selected.  The  difference  has  a negligible 
effect  on  the  wall  temperature  history. 
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APPENDIX  B 


Heat  and  Mass  Tranafer  Express ions 


1«  Turbulent  Natural  Convection 


Lav** 


Using  the 
equations  8.6Q| 
unit  width  (ifl 


j^Vlts  r^ve^  by  Rohsenow  and^^^ 
8. si),  it  can  be  shown  that^^^ 
Win  *-h*  turbulent,  free  connect io^ 
at  *11)  section  at  a height  t above  t 


(1961,  page  204, 
kMss  flow  rate  per 
^L'dary  layer  in  anv 
b^BM  is  given  by 


general  f 


Crashoff  nuaber 


' * Prandtl  nunber  of  the  fluid 


iV«  the  aubscript  • refers  to  the  condition  far  fr« 


Table  A-1 


Results  of  the  Tank  Heating  %d.th  Stratification 
of  Liquid  Considered 


t Hall 

Tine  Temperature 

(Secs)  (C) 


T 

s 

Liquid 

Surface 

Temperature 

(C) 


S 

Depth  of 
Stratified 
Layer 
(ft) 


Physically  unrealistic 
predictions  because 
boiling  has  not  been 
included  and  the  wall 
temperature  exceeds 
the  boiling  point  of 
most  cargoes. 


-01 

1 ,667b 

01 

l,667b 

01 

3 F.-Ol 

1,282F 

01 

3, Jibt 

01 

2,391b 

01 

1 ,54bl.-e2 

?,563e 

01 

4.90lt 

01 

2,843F 

01 

4,8/7( -82 

3,84»>f 

01 

6,3J6£ 

01 

3,Jb0E 

01 

e,340t-H2 

3,aebE 

01 

6,379t 

01 

3,3b9E 

01 

8,4691 -/2 

3.92bE 

01 

6,422t 

01 

3,374b 

01 

8,5941 -02 

3,96bE 

01 

6,464b 

01 

3,388b 

01 

8,7201-0? 

d.OHSt 

01 

6,ft07t 

01 

3,404E 

01 

8,846e-02 

4,H4St 

01 

6,b49b 

01 

3,418b 

01 

8,9731 -02 

4,88b( 

01 

6,691e 

01 

3,433E 

01 

9,098t-02 

4,165P 

0] 

6,6/4^ 

01 

3,463F 

01 

9,3531-02 

4,245fc 

01 

6,7b7b 

01 

3,493b 

01 

9,6101-02 

4,32bE 

01 

6,848F 

01 

3,b22b 

01 

9,8671 -02 

4,4tfl6E 

01 

h,92U 

01 

3,552b 

01 

1 ,013F-01 

4.56br 

01 

7,084b 

01 

3,610E 

01 

1 ,M65t-01 

4,726fc 

01 

7,2441 

01 

3,068b 

01 

1,1181-31 

4,886L 

01 

7,40Pe 

01 

3,726F 

01 

1 , 1 7 1 1 -(.*  1 

5,(tf46t 

01 

7,S57b 

01 

3,783E 

01 

1 ,2251-01 

5, Jb7f 

01 

7,8b2b 

01 

3,89bb 

01 

1,3361-01 

b,687t 

01 

8,lb9E 

01 

4,006b 

01 

l,447F-,n 

b.dMet 

01 

8, 447b 

01 

4,ltbb 

01 

1,5611-01 

6,328e 

01 

e,727E 

01 

4,222b 

01 

l,677f-01 

b,969fc 

01 

9.26?b 

01 

4,431E 

01 

1,9131-01 

7,b|kJt 

01 

9,/6bE 

01 

4,b32b 

01 

2,1551-01 

01 

1 ,024t 

02 

4,82bE 

01 

?,4|tt3t-01 

8,b91k 

01 

1 ,068t 

0? 

5,0llt 

01 

2,6561-01 

J ,4’17f 

02 

M49t 

02 

5,362e 

01 

3, I74t-0l 

l.ldbf 

02 

1,2191 

0? 

5,686E 

01 

3,7061-01 

l,274fc 

02 

1 ,2H0t 

02 

b,98bt 

01 

4,2501  -in 

l,«'d2£ 

02 

1,333b 

02 

6, 261b 

01 

4,81 2t-toJ 

1 ,bh8£ 

02 

l,419b 

02 

6,7b0F 

01 

5,9271-01 

1 ,9l4t 

02 

1 ,483b 

02 

/,168f 

01 

7, 0681 -01 

TTTTTr 

02 

r.bilb 

02 

/,b2/b 

01 

8,2|9l.0n 

?,d27t 

02 

1 ,b66t 

02 

7,836f 

01 

9,3731-01 

2, 6831 

02 

1 ,b93t 

02 

H,10bF. 

01 

1,0531  01' 

2,9421. 

02 

1 ,6131 

02 

b,341F. 

01 

1,1681  01 

3,l9bt 

02 

1 ,627f 

02 

6,549F 

01 

l,2e2.t  0»i 

3,4b?£ 

02 

1,639b 

02 

8,73bF 

01 

1,3961  00 

3,7H9t 

02 

1 ,647f. 

02 

8,902E 

01 

1,5091  00 

3,96bt 

02 

1,654b 

02 

9,052F. 

01 

1,621F  01' 

4,221F 

02 

1 ,660b 

02 

9, 190r 

HI 

1,7331  (“0 

4,4/81 

02 

1 ,6b4t 

02 

9,317b 

01 

1,8431  01- 

2 

Incident  - heat  flux  « 35,000  Btu/hr/ft 

T^  ■ wall  temperature  in  contact  with  the  liquid 
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APPENDIX  B 


Heat  and  Mass  Transfer  Expressions 


Turbulent  Natural  Convection 


Turbulant  - Natural 
Convection  Boundary  Layer 


FIGURE  B-1  SCHEMATIC  REPRESENTATION  OF  THE  FLOW 
IN  A NATURAL  CONVECTION  HEATING  OF  A 
FLUID  BY  A VERTICAL  WALL 


using  the  results  given  by  Rohsenow  and  Choi  (1961,  page  204, 
equations  8.60  and  8.61),  it  can  be  shown  that  the  mass  flow  rate  per 
unit  width  (M  ) in  the  turbulent,  free  convection  boundary  layer  in  a 
general  fluid  at  any  section  at  a height  z above  the  base  is  given  by 


M'(z)  - 0.1  (p  v) 


Pr®^^^  tl  + 0.494  Pr^^^] 


where  Gr  ■ Grashoff  number 
z 


- TJz- 


and  Pr  ■ Prandtl  number  of  the  fluid 

where  the  subscript  » refers  to  the  condition  far  from  the  wall. 


Similarly,  from  equation  8.62b  of  Rohsenow,  we  have,  for  the  average 
Nusselt  nun^er  for  heat  transfer  over  a wall  height  of  z. 


0.02A6 


Gr"-^  Pr’/^5 

z 


2/3 

[1  + 0.494  Pr  ' ] 


(B-3) 


Therefore,  the  average  heat  flux  from  the  wall  to  the  fluid  over  a 
height  z is 

T-  - K(T  - Tj 

ql’  - h(T  - T ) » Nu  ^ 

X W “ z z 

We  also  define  a quantity  called  the  characteristic  mean  heat  flux. 


0.4 


K T 

- 0.0246 


get  H- 


Pr 


7/15 


2/3 

[1  + 0.494  Pr  ] 


0.4 


(B-5) 


Physically,  represents  that  average  heat  flux  which  would  result 

over  a height  H when  the  wall  temperature  Is  twice  the  absolute  tempera- 
ture of  the  bulk  of  the  fluid. 

The  quantity  T defined  in  Appendix  A,  equation  A-13,  is  therefore 
the  ratio  of  the  above  characteristic  mean  heat  flux  q|]j^  over  the  liquid 
height  to  the  heat  flux  from  fire. 

2.  Radiation  View  Factor  Between  Two  Normal.  Rectangular  Surfaces 


FIGURE  B-2  CONFIGURATION  FOR  RADIATION  VIEW  FACTOR 


B-2 


With 


. _v 
h - „ 

b 

W 

the  view  factor  is  given  by 


. i_  fh  tan"^  0+  b tan~^(^)-  + b^  tan'^ 

irh  1 + 

1 r,  2 h^(l  + + b^)  2 b^d  + + b^) 

(1  + h^)(h^  + b^)  (1  + b^)(h^  + b^) 

^ SI*  Aa-±|^i1 
(1  + h + b ) J 


(B-6a) 


(B-6b) 


(B-7) 


The  above  formula  is  taken  from  Wiebelt  (1966),  page  259. 

The  total  net  heat  transfer  rate  from  surface  1 to  surface  2 is 
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Derivation  of  Energy  Conservation  Equation  for 
the  Gas  In  the  Ullage  Volume 


Objective 

The  objective  of  the  analysis  presented  In  this  appendix  Is  to 
obtain  expression  for  the  change  with  time  of  the  temperature  of  the 
gas  In  the  ullage  volume  of  the  tank  when  exposed  to  fire  and  subjected  to 
mass  Input  and  venting. 

Analysis 

In  order  to  obtain  the  gas  temperature  variation,  we  write  the 
energy  conservation  equation  subject  to  the  following  conditions: 

• The  kinetic  energy  of  flow  is  neglected  in  comparison  with 
the  thermal  flow  during  venting. 

• Gases  behave  as  perfect  gases  which  are  mechanically  mixed  (no 
chemical  reaction  takes  place  between  air  and  vapor). 

• Direct  radiative  heating  of  gas  due  to  thermal  radiation  from 
the  wall  Is  small  compared  to  the  gas  heating  by  natural  con- 
vection. 

• The  gas  Is  In  a perfectly  mixed  state  always;  that  Is,  the  gas 
stratification.  If  any,  is  neglected. 

Figure  C-1  is  a schematic  representation  of  the  thermodynamic  system 
(of  constant  volume)  containing  the  ullage  gas. 


Total  Internal  Energy  > E 

Me  6^ 

FIGURE  C-1  SCHEMATIC  REPRESENTATION  OF  THE  THERMODYNAMIC  SYSTEM 
CONTAINING  THE  ULLAGE  GAS 


i 


L 


c-1 


with 


e ■ specific  internal  energy  of  individual  gas  specie  (l.e.,  per 
unit  mass) 

h « specific  enthalpy  of  the  individual  gas  specie 

E » total  internal  energy  of  the  gas  in  the  ullage  volume  at 
any  Instant 

We  write  the  energy  equation  to  this  open  system  (see  Keenan  (1970) , 
page  32]. 

^ « Q + . - M h 

dt  c i 1 o o 

where  h and  h are  the  gas  enthalpies  at  inlet  and  outlet  per  unit 
i o 

mass  of  the  gas  mixture. 

E - M[Y  e + Y e ] 
a a V v 


(C-1) 


(C-2) 


where  Y is  the  mass  fraction  of  air  in  the  ullage  volume  and  Y is 
a V 

the  mass  fraction  of  vapor.  Also, 

Y + Y - 1 
a V 


(C-3) 


Differentiating  Equation  (C-2)  and  substituting  the  results  in  (C-1),  we  get 
(in view  of  C-3). 

• • * • ■ • — • • _ 

M[Y  e +Ye]+MY[e  - e]+M[Ye  + Ye]-Q  +Mh  - Mh  (u-4) 
‘aa  vv  a a v aa  vv  ^c  ii  oo  ' 


vdiere  the  "dot"  over  the  symbol  represents  its  rate  of  change  with  time. 

From  the  mass  balance  equation,  we  have 

• • • 

M - M - M 
1 o 


(C-5) 


Substituting  C-5  in  C-4  and  rearranging  after  substituting  for  e its 
V dT 

equivalent  c we  have 

at 

V V dT  * * — • — 

M[Yc  + Yc]^-Q  +Mrh-Ye  - Ye]-M[h  - Ye  - Ye] 
a a V V dt  c 11  a a v v o o a a v v 


(C-6) 


For  Y^,  we  substitute  from  equation  10  in  equation  C-6  and  get 


C-2 


M[Y  + Y c^]  4r  * Q + M Jh  - Y e -Ye  + Y e -Ye] 

aa  vvdt  c il  aa  vv  aa  av 


-M[h  - Ye  - Ye] 
o o a a V V 


Noting  that  in  the  inlet  gas  there  is  only  vapor  and  that  outlet 
gas  contains  both  air  and  vapor. 


h j , * h , 

V,  inlet  saturated 


h - Y h (T)  + Y h (T) 
o a a V V 


we  have 


M[Y  c''  + Y c ] 4^  = Q + M fh  - h (T)  + ^]  - M [Y  (h  - e ) 
a a V V dt  c i sat  v p o a a a 


+ Y (h  - e )] 

V V V 


1.  e. 


l.e. 

• • 

-R  -| 

dT  . 

Qc  + «i 

_li  T + c^(T.  - T) 

U V D 

L V 

- M R T 
o u 

^ + 1 
P P 

La  V J 

dt 


V V 

MtY  c + Y c ] 
a a V V 


It  is  also  recalled  that 
,P 


c' 

c 


. p V 
and  c*^  - c 


V/  IS  u 

c (Y  - 1)  - — 


(C-7) 


(C-8a) 

(C-8b) 


V .fRT  I-Ty  YT 

M[Y  c + Y c'^]  • Q + M.  1 — + c^’d,  - T)  - M — + -^  R 1 

a a V V dt  c v b J o J u 


(C-9) 


(C-10) 

(r-11) 


Equations  C-10  and  C-11  apply  equally  to  air  and  vapor  (and  have  to  be 
used  with  appropriate  subscripts). 


It  is  also  noted  that  the  mean  values  of  specific  heat  (mass)  and 
molecular  weights  in  a mixture  of  gases  (say,  containing  air  and  vapor) 
are  expressed  by 
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APPENDIX  D 

Non-Dlmenslonalization  and  Solution  Procedure 
for  the  Governing  Equations 

In  this  appendix,  the  equations  derived  In  the  analysis  of  tank 
heating  and  pressurization  are  written  In  dimensionless  form  so  as  to 
be  amenable  to  numerical  solution  on  a computer.  Also  indicated  briefly 
is  the  method  of  solution. 


We  first  define  the  following  non-dlmenslonallzing  parameters, 
WHl 


M 


ch 


characteristic  mass  generation  rate 


(D-la) 


represents  the  vapor  generation  rate  that  would  occur  If  all  of 
the  heat  incident  on  the  liquid  wetted  wall  face  is  absorbed  by  the 
liquid  during  the  boiling  process. 


M 


Pin" 
In  R 

T 

- in 
‘^In 


initial  mass  of  the  vapor-air  mixture  in 
the  ullage  volume 


where 


Y Y 

1 _ a, in  ^ v,ln 

— u u 

Pm 


L p c T_, 

^ _ w w w in 

ch 


characteristic  time 


Physically,  t^j^  represents  the  time  in  which  the  wall  temperature  would 
increase  by  its  initial  temperature  if  all  of  the  incident  radiation 
is  absorbed  by  the  wall  itself. 

o 

o’  ■ — ■ dimensionless  radiative  flux  from  the  wall 

q^  initially 

O.A 


Nu  - 0.0246 


1 


8 PT,  H 


3n 


In  V 


(Pr.) 


7/15 


II  + 0.494 


0.4 


(D-lb) 


(D-lc) 


(D-ld) 


(D-le) 


(D-lf) 


I>-1 


Nu*  Is  a characteristic  mean  Nusselt  number  for  turbulent  convective 
heat  transfer  In  air  over  a height  when  the  temperature  difference 
between  the  wall  and  air  Is  equal  to  the  Initial  temperature.  The 
product  of  Equation  (D-lf)  Is  unity  for  a perfect  gas. 

p*  - 

n . _V  1: 

H-  T.  (R  /u  ) 

X in  u V 


C A 1/2  p p 
d ^^In  ^In 


Physically,  A represents  the  ratio  of  the  outflow  rate  through  the  vent 
when  there  Is  a pressure  drop  across  the  vent  equal  to  the  Initial 
pressure  In  the  tank  to  the  characteristic  evaporation  rate  defined 
in  equation  D-la. 


e * T /T  • liquid  boiling  temperature 
b b^  in  Initial  vapor  temperature 

e ■ T /T  « 1 « bulk  liquid  temperature 
D b'  In  Initial  temperature 


r 


c 


ch 


^ ^amblent^**ln 

heat  flux  that  can  be  transferred  in  air  by  conduction 
C ■ — ; — ■ for  a temperature  difference  of  T^jj  over  a length  H 

{{  q"  — i 

V f Incident  fire  heat  flux 


Non-dimensional  Variables 

T * t/t  . ~ dimensionless  time 
ch 

6 ■ T/T.  ■ dimensionless  vapor  temperature 

in 

e a T /T.  ■ dimensionless  wall  temperature 

w w in 

_ w/%.  _ mass  of  vapor-alr  In  the  ullage  volume 

m ■ n/  n . " \ •*— — — — 

In  mass  of  gases  initially 


**ln 


(D-lg) 

(D-lh) 

(D-11) 

(D-lj) 

(D-lk) 

(15-11) 

(D-lm) 

(D-ln) 

(D-lo) 

(D-2a) 

(D-2b) 

(D-2c) 

(D-2d) 


dimensionless  pressure 


(t^2e) 


“i 

— - dimensionless  vapor  mass  inflow  rate  by  liquid 

M . boiling 
cn 

M 

■ dimensionless  mass  outflow  by  venting 


M 


ch 


q ~ q"  /q"  ■ dimensionless  natural  convection  heat  flux 
c c f 

q_  . ■ q"  /q"  ■ dimensionless  radiative  flux  leaving 
Rad  rad  f -i  i 

the  wall 


Using  the  above  definitions,  equations  1 through  13  are  written  in 
dimensionless  form  as  shown  below: 

Equation  2 becomes 

q - o'  e ^ 

^Rad  w 


Equation  3 becomes 


q - C Nu*  (0  - 0) 

^c  w 


7/5 


(D-2f) 

(D-2g) 

(D-2h) 

(D-2i) 


(n-3) 


(D-4) 


The  energy  equation, 
de  - » 

di  p* 


Equation  (13),  becomes 

(i!  + .j{e  + - «))  - "„(V*  * '’v>  »1 

bTy  + c Y ] ^ (D-10) 

a V 


Finally,  we  have 

Y + Y - 1 
a V 


Initial  Conditions 


a,  in 

e - e 

w 


/ 

at  T 

1 


0 


(D-11) 


(r'-i2) 


Method  of  Solution  o'  tlic  E aations 

There  are  essentially  four  unknowns;  namely,  6 , Y and  0, 

W fl 

and  four  equations  describing  their  time  rate  of  change;  namely,  D-5, 

D-8,  D-9,  and  D-10  respectively.  However,  these  differential  equa- 
tions are  coupled  and  nonlinear.  Only  a numerical  solution  can  be 
obtained. 

The  method  used  is  the  Hamming  Predictor  Corrector  Method  enunciated 
by  Hananlng (1973) . The  program  used  is  the  one  available  in  the  Scientific 
Systems  Subroutine  Package  developed  by  IBM. 


D-4 


APPENDIX  E I 

Analytical  Analysis  of  the  Pressure  Rise  Problem  | 

A Simple  Model j 

I 

Objective 

The  objective  of  the  analysis  developed  in  this  appendix  is  to 
obtain  a closed  form  solution  to  the  variation  of  tank  pressure  as  a 
function  of  time  using  certain  simplifying  assumptions. 

Introduction 

The  calculation  of  pressure  rise  with  time  resulting  due  to  the 
tank  being  exposed  to  a fire  Is  a complicated  problem.  For  most  con- 
ditions encountered  In  practice,  the  solution  of  the  equations  can  be 
obtained  only  by  numerical  techniques  using  a computer. 

If  one  analyzes  the  effect  of  various  physical  phenomena  that 
cause  the  pressure  rise,  one  finds  that  the  main  effect  is  due  to 
the  tremendous  vapor  Inflow  Into  the  ullage  space  due  to  the  liquid 
boiling.  In  fact,  most  of  the  pressure  rise  in  the  Initial  stages  Is 
only  due  to  this  mass  Influx.  The  thermal  effect  (l.e.,  the  effect 
of  Increase  In  temperature  of  the  gas  on  the  gas  pressure)  comes  Into 
play  at  a much  later  stage. 

Therefore,  In  light  of  the  above  observation,  the  analysis  can  be 
simplified  to  a great  extent  resulting  In  closed  form  solution  of  pres- 
sure as  a function  of  time.  In  effect,  then,  the  analysis  will  be  an 
Isothermal  analysis. 

The  following  assumptions  are  made  In  this  analysis: 

• All  gases  are  perfect  gases. 

• The  temperature  of  the  gases  remains  at  a constant  (Initial)  value. 

• The  mass  inflow  rate  of  vapor  into  the  ullage  volume  by  boiling 
of  liquid  la  constant. 

0 The  vapor  has  the  same  snlecular  characteristics  as  that  of  air. 

s Volume  of  vapor  space  is  a constant. 

0 Initial  pressure  is  equal  to  the  aad>ient  pressure. 
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Analysis 


Let  * constant  Inflow  rate  of  vapor  Into  the  ullage  space 
We  now  write  the  mass  balance  equation  on  the  ullape  vapor: 


f-  «i  - - ■’a.*’'’ 

t ( I ' 

Inflow  venting 


with 


M * * mass  of  gas  in  the  ullage  space 

Rx 


M 


Equations  E-1,  E-2,  and  E-3  have  to  be  solved  for  p (the  pressure) 
as  a function  of  time  with  the  Initial  condition 

P - Pamb  ^ ° 

In  order  to  obtain  closed  form  solution,  we  define  the  following 
parameters: 


amb 


In  RT 


In 


M,  - P,  V 
In  In 


t - 

”l 


[f  ■ ^ because  T is  constant 

Pamb  ”ln 


Using  the  parameters  defined  In  equations  E-5a  through  E-5f,  we 
rewrite  equation  E-1  as 

if-  

dT 


(E-1) 

(E-2) 

(K-3) 

(E-4) 

(E-5a) 

(E-5b) 

(E-5c) 

(E-5d) 

(E-5e) 

(E-5f) 


1 - A 1) 


(E-6) 


The  exact  solution  to  Equation  (E-6)  can  be  obtained  using  standard  re- 
sults for  Integrals  given  In  Gradshteyn  and  Ryzhlk  (1973),  page  107, 
equation  2.441-3.  The  solution  to  Equation  (E-6)  with 


4 ■ cosh  (E-9) 


N(»ffiNCLATURE 


a 

A 


b 

B 


c 


D 

E 


F 

Gr 

8 

«L 

H 

V 

h 


K 


! 


> Coefficient  In  Equation  (22) 

■ Cross  sectional  flow  area  of  the  vent 

> Liquid  surface  area  » BW 

■ Plating  width 

m Beam  width  of  the  tank 

> Coefficient  of  discharge  of  the  vent 

> Mass  specific  heat  at  constant  pressure 

> Mass  specific  heat  at  constant  volume 

- Specific  heat  of  liquid 

- Vapor  to  air  specific  heat  ratio  (Equation  (D-lm)) 

* Diffusion  coefficient 

~ Internal  energy  content  of  the  gases  in  the 
ullage  volume 

Also  used  for  Modulus  of  Elasticity  in 
Chapter  II 

> Specific  internal  energy  of  individual  gas  specie 

- Radiation  view  factor  (Equation  (B-7)) 

> Grashoff  number  (Equation  (B-2)) 

* Acceleration  due  to  gravity 

■ Depth  of  liquid  above  the  water  line 

■ Depth  of  vapor  space 

> Turbulent,  natural  convection  transfer  coeffi- 
cient 

Also  used  for  specific  enthalpy  in  Appendix  C 
(Equation  (C-1)) 

Also  used  for  plating  thickness  in  Chapter  II 
" Thermal  conductivity  of  air 


(hr-ft^/lb)^"^  psl 

ft^ 

ft^ 

ft 

ft 

Btu/lbm/*F 

Btu/lbm/'’F 

Btu/lbm/'F 

Btu/lbm/®F 

ft^/hr 

Btu 

psl 

Btu/lbm 


2 

ft/sec 

ft 

ft 

Btu/hr/ft^/*F 

Btu/lbm 

ft 

Btu/hr/ft/*F 
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L ■*  Thickness  of  steel  wall 
w 

M » Mass  of  air  and  vapor  In  the  ullage  volume 


M •*  Mass  of  air  In  the  ullage  volume 
d 


« Mass  rate  of  liberation  of  vapor  due  to  liquid 
boiling 


M 


evap 


Mass  rate  of  liquid  evaporation  at  the  surface 
due  to  direct  radiation  from  the  hot  wall  In 
contact  with  the  vapor 


M. 


~ Vapor  mass  Inflow  rate  Into  ullage  space  due  to 
liquid  boiling  and  evaporation  at  the  surface 


M.  * Initial  mass  of  vapor  in  the  ullage  volure 
In 

M ” Mass  outflow  rate  through  the  vent 
o 

M * Mass  of  vapor  In  the  ullage  volume 

hH 

V 

Nu  “ Musselt  number  “ — 

Nu*  ” Characteristic  Nusselt  number  (see  equation  D If ) 


Pr  “ Prandtl  number  of  air  or  vapor 

(p  ■ Dimensionless  tank  pressure  ■ (equation  D-2e) 

p > Pressure  in  the  tank 

p^  ■ Partial  pressure  of  air  in  ullage  volume 
Pamb  “ Ambient  pressure 

p ~ Initial  pressure  of  gases  In  the  tank 

p > Steady  state  tank  pressure 

68 

p^  * Partial  pressure  of  vapor  In  the  ullage  volume 

p"  ■ Pajni,/Pin  (Equation  (D-ln)) 

P > Gauge  pressure  at  the  tank  failure  limit 

ttUX 

Q " Total  heat  Input  rate  Into  the  gas  by  natural 
convection 

q - Dlmenslonleaa  liquid  convection  heat  flux 
(Equation  (A-'7e}) 


ft 

Ibm 

Ibm 

Ibm/hr 


Ibm/hr 

Ibm 

Ibm/hr 

Ibm 


psl 

psl 

psl 

psl 

psl 

psl 

pslg 

Btu/hr 
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• 2 
> Mean  heat  flux  from  the  wall  to  the  fluid  in  Btu/hr/ft 

contact,  due  to  turbulent  free  convection 

(equation  B-4) 

q"  ■ Mean  turbulent  heat  flux  for  natural  convec-  Btu/hr/ft 

tlon  in  liquid 


q'^  * Heat  flux  at  the  outer  surface  wall  due  to  the 

Incoming  radiation  from  a fire 

q^^j  ■ Heat  flux  leaving  the  wall  (in  contact  with  the 
vapor)  by  radiation 

R > Individual  gas  constant 

R <■  Universal  gas  constant 

u 

S > Depth  of  the  stratified  liquid  layer  at  any  time 

T - Temperature  of  the  gases  In  the  tank 

T,  ■ Bolling  temperature  of  liquid 
b 

“ Surface  temperature  of  liquid 
” Wall  temperature 
t » Time 

t . • Characteristic  time  (see  equations  A-7b  and  D-ld) 

cn 

V ■ Volume  of  the  ullage  space 

W - Width  of  tank  exposed  to  fire  radiation 

Y - Mass  fraction  of  the  air  or  vapor  specie  in 

the  ullage  gas 

z - Vertical  coordinate 


Btu/hr/ft^ 

Btu/hr/ft^ 


Btu/lbrn/'F 

I 

Btu/lbmole/“R 

ft 

“R 

*R 

•r 

•r 

sec 

sec 

ft^ 

ft 


ft 
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GREEK  LETTERS 


m Volumetric  expansion  coefficient  of  a fluid 
(for  perfect  gas  6 ■ 1/T) 

cP 

- Specific  heat  ratio  of  a gas  - — 

c’’ 

« Parameter  defined  In  equation  E-5f. 

- Parameter  defined  In  equation  A-13 

■ Dimensionless  liquid  stratification  de^tH.  Also 
a parameter  defined  In  equation  D-lo 

“ Dimensionless  temperature  ■ 

■>  Heat  of  vaporization  of  liquid 

■ Parameter  defined  In  equation  A-7h 

■ Molecular  weight  of  a gas  specie (subscripted) 

■ Mean  molecular  weight  of  air  vapor  mixture 
(see  equation  C-12) 

% 

■ — “ Ratio  of  vapor  to  air  molecular  weight 

a (equation  D-lg) 

« Parameter  defined  In  Equation  (18) 

> Kinematic  viscosity  of  the  fluid 

■ Density  of  fluid  (subscripted)  or  wall 
* Stefan  Boltzmann  Constant 

Also  used  for  yield  strength  In  Chapter  II 


* Dimensionless  radiative  heat  flux 


T ■ T — “ dimensionless  time  (see  equation  A"7c) 

ch 

♦ - Pressure  parameter  defined  in  equation  E-9 

u “ Parameter  defined  In  equation  A-7f 

fi  ■ Parameter  defined  i’  D-lh 


Btu/lbm 


ft^/hr 


Ibm/ft 


Btu/hr/ft^/*R^ 
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SUBSCRIPTS 
a 

B ■ Bulk  liquid 

b * Boiling  condition 

c « Convection 

ch  * Characteristic  value 

e ■ Evaporation 

f * Flame 

i » Inlet  (ullage  volume)  condition 

in  * Initial  condition 

L = Liquid  condition 

o » Outlet  (vent)  condition 

s ■ Surface  of  liquid 

ss  * Steady  state 

sat  * Saturation 

V « vapor 

® • for  conditions  far  from  the  wall 

SUPERSCRIPTS 

. ■ Time  rate  of  change 

' * Per  unit  length 

" ■ Per  unit  area 

p « At  constant  pressure 

V - At  constant  volume 
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